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We present a large-scale molecular dynamics study of nematic-paranematic interfaces under shear. 
We use a model of soft repulsive ellipsoidal particles with well-known equilibrium properties, and 
consider interfaces which are oriented normal to the direction of the shear gradient (common stress 
case). The director at the interface is oriented parallel to the interface (planar). A fixed average 
shear rate is imposed with Lees-Edwards boundary conditions, and the heat is dissipated with a 
profile-unbiased thermostat. First we study the properties of the interface at one particular shear 
rate in detail. The local interfacial profiles and the capillary wave fluctuations of the interfaces are 
calculated and compared with those of the corresponding equilibrium interface. Under shear, the 
interfacial width broadens and the capillary wave amplitudes at large wavelengths increase. The 
strain is distributed inhomogeneously in the system (shear banding), the local shear rate in the 
nematic region being distinctly higher than in the paranematic region. Surprisingly, we also observe 
(symmetry breaking) flow in the vorticity direction, with opposite direction in the nematic and the 
paranematic state. Finally, we investigate the stability of the interface for other shear rates and 
construct a nonequilibrium phase diagram. 



I. INTRODUCTION 

Complex fluids often show peculiar behavior under 
sheari Particularly remarkable phenomena are observed 
in fluids with anisotropic components, which can exhibit 
local or global orientational orderiS*^ Shear flow tends 
to orient the materials, which usually reduces the shear 
viscosity (shear thinning) ^i^i^ but can also have the op- 
posite effect In the vicinity of equilibrium phase tran- 
sitions, shear flow may shift the boundaries and enforce 
ordering transitionsAiSiiiiiSiiiii It may even induce the 
formation of new structures that are inhomogeneous in 
space or time, such as multilamellar vesicle phases in 
surfactant systemsi^ti^ or tumbling states in liquid crys- 
talline polymers i2ii2*i2iSai2i 

In recent years special attention has been given to sit- 
uations where phases "coexist" under shear. The inter- 
est was originally raised by the "spurt effect" in poly- 
mer melts, i. e., the observation that the flow rate in 
a pipe may change discontinuously as a function of the 
applied pressure differencei^^ This was attributed to a 
nonequilibrium phase transition from an entangled poly- 
mer phase with high viscosity to a highly aligned phase 
with low viscosity The instability results from 
a nonmonotonicity in the (hypothetical) homogeneous 
stress-strain flow curve: In a homogeneous system, the 
stress would first grow with growing strain rate and then 
decrease beyond a certain critical value. Since such a 
flow curve is mechanically unstable, the system becomes 
inhomogeneous. The stress plateau in the stress-strain 
flow curve was interpreted as the signature of a region 
where phases with different shear rates coexist. The sep- 
aration of a fluid into two states with distinctly different 
shear rates is commonly called shear banding. 

Meanwhile, shear banding has been observed and 
even directly visualized in various systems, most notably 
wormlike micellesSSiSSiSLSiSS and rodlikc virusesiSi More 



generally, two qualitatively different types of phase sep- 
aration are possible in shear flowi^ These are illustrated 
in Fig. n for the Couette geometry. In the shear banding 
case described above, the fluid separates in the direction 
of the flow gradient, the phases are subject to the same 
stress, and the shear rates differ (common stress case). 
Alternatively, the fluid can also phase separate in the di- 
rection perpendicular to both flow and flow gradient, the 
vorticity direction. The two phases then share the same 
shear profile, but the stress differs (common strain case). 

The phase separation mechanism described above is 
driven by mechanical (hydrodynamical) forces. However, 
phase separation and shear banding can also occur as 
a result of a thermodynamic phase transition - either a 
real transition that will also occur in the absence of shear 
flow, or a transition between phases that are metastable 
("hidden") in equilibrium. Shear banding is often associ- 
ated with ordering phenomena and sometimes occurs in 
the vicinity of a thermodynamic phase boundary. Ther- 
modynamic forces then contribute to the nonmonotonic 
stress-strain relationship which gives rise to the mechan- 
ical instability. Thus the mechanical and the thermo- 
dynamical phase separation cannot always be separated 
unambiguously 1^ In fact, theoretical treatments of shear 
banding are often based on free energy expressions for 
phase separating systemsi^2iSi2li^ 

One should, however, note an important difference be- 
tween mechanically and thermodynamically driven phase 
separation^ The signature of mechanically driven phase 
coexistence is the existence of a plateau in the stress- 
strain flow curve. In the common stress case, the plateau 
is horizontal, in the common shear case, it is verticali^ If 
one increases the average shear rate/average stress in the 
plateau region, i. e., the coexistence region, the width 
of the coexisting bands readjusts, but the local stress 
and the shear rate in each phase remain the same. For 
thermodynamically driven phase separation, however, a 
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free readjustment is only possible if solely nonconserved 
quantities are involved in the phase transition. If the two 
coexisting phases differ with respect to the density of a 
conseraerf quantity (e. g., a concentration) , the overall av- 
erage of this density imposes an additional constraint. As 
a result, both the band width and the coexisting phases 
must be readjusted, and the "plateau" acquires a finite 
slopeMi^i 

This raises an important question: How does the sys- 
tem determine the location of the plateau and select the 
coexisting states? According to Olmsted and coworkers, 
the selection criterion is the requirement that a steady 
interface between stable coexisting homogeneous states 
exists Hence the understanding of the interface 
provides the key to the understanding of the properties 
of the material under shear. 

In the present work, we present a molecular dynam- 
ics simulation of a nematic/isotropic (more precisely, 
paranematic) interface under shear, in a situation where 
the interface is already present in the absence of shear. 
The phase separation is thus driven by thermodynamic 
forces. The equilibrium properties of isotropic/nematic 
interfaces have been studied by computer simulations in 
a number of systemsi^2i^2i^2*^ii^2i^2ili Our study is de- 
signed such that we can compare the equilibrium and 
the nonequilibrium structure of the interface in as much 
detail as possible. Therefore, we choose a model of 
soft ellipsoids, which we have investigated extensively in 
the pastpi'^2iiSiiSiiLi£ and used for large scale simula- 
tions of equilibrium interfaces4i*^ These interfaces are 
now subjected to shear. Previous particle-based simu- 
lations of liquid crystals under shear have focussed on 
the structure and the response functions in homogeneous 
phasesi'?i^?i^'??'^^i^^n'^?i^'^^i-'^^i'?^''?Ti'?? To our best knowledge, 
our simulation is the first particle-based nonequilibrium 
simulation of a liquid crystal interface under flow. 

Our paper is organized as follows: In the next sec- 
tion, we introduce the model and discuss the simulation 
method. The results are presented in section lllll We first 
discuss the stability of the interface with increasing shear 
rate, and then consider in detail the interfacial structures 
at one particular shear rate. We summarize and conclude 
in section Hvl 
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II. MODEL AND METHOD 

We consider a system of soft repulsive ellipsoids with 
length to width ratio ai/cro = 15. Every ellipsoid is char- 
acterized by the position of the center of mass r and a 
unit vector u pointing in the direction of the long axis. 
The interaction potential between two ellipsoids i and j 
separated by the center-center vector is given by 
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where Sij = (rij — cr(fy, u^, Uj) -I- (To)/o'o is a scaled and 
shifted distance and a approximates the contact distance 
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FIG. 1: Possible orientations of an interface: common strain 
(left) and common stress (right) geometries. Also indicated 
is the reference system used throughout this paper. In this 
paper, only the common stress case is investigated. 



of two ellipsoidsr^ 
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The particles have the mass m and the moment of inertia 
/ = 50 mcTQ. For convenience, we choose the units such 
that ao ~ eo = m = 1. This defines the time unit r = 
(To i/m/eo. The temperature was fcsT = eg. 

In our previous work^ we have performed an equilib- 
rium molecular dynamics simulation of a system with 
115200 particles in an elongated box with side ratios 
(1:2:1) and periodic boundary conditions, at a density in 
the coexistence region, p = O.OIT/ctq. The system thus 
contained an isotropic and a nematic slab, separated by 
two interfaces. The director in the nematic phase was 
found to align parallel to the interface (planar). Within 
that plane, the direction of the director is of course ar- 
bitrary, but the initial equilibration stage was conducted 
such that it pointed roughly in the direction of one box 
side. 

The final configurations from these simulations were 
used as starting configurations in the present work. Shear 
flow was enforced by applying Lees-Edwards boundary 
conditions r^iSiiS i. e., the periodic boundaries in one di- 
rection move at constant speed. This minimizes surface 
effects, but at the expense of constantly pumping energy 
into the system. In order to dissipate the latter, we have 
coupled the molecular velocities to a Nose-Hoover ther- 
mostatfS taking care that it acts only on the thermal 
fluctuations of the velocities, i. e., the velocities minus 
the streaming velocity of the fluid. The streaming ve- 
locity profile was determined from the simulation. In a 
nonuniform system like ours, using such a profile unbi- 
ased thermostaliS is obviously important. However, we 
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have noted in preliminary studies that we get practically 
the same results with a profile biased thermostat, where 
the actual velocity profile is approximated by a constant 
gradient profile. Hence the results do not seem to depend 
sensitively on the details of the thermostat. 

In the present paper, we discuss interfaces in the com- 
mon stress geometry (Fig. ^ left). Two situations were 
considered: In the first, termed "flow-aligned" case, the 
initial configuration was turned such that the director in 
the nematic phase points in the direction of the flow, 
i. e., the x direction in Fig. ^ In the second, "log- 
rolling" case, the initial configuration was turned such 
that the director points in the direction of the vorticity, 
i. e., the z direction. We have first examined the sta- 
bility of the interface in the flow aligned case for global 
strain rates up to 7 = Avx/Ly ~ O.l/r. (Au^ is the 
streaming velocity difference between both sides of the 
box, as imposed by the Lecs-Edwards boundary condi- 
tions, and Ly the box length in the y direction.) To 
this end, we have carried out a scries of simulations of 
smaller systems with 7200 particles in a box with side 
ratios (1:8:1). Then we have focussed on the strain rate 
7 = 0.001/t, and performed extensive simulations of in- 
terfaces with initial flow-aligned and log-rolling setup. 
In both cases the system was first "equilibrated" over 
2 million MD steps, and data were then collected over 
the next 2 million MD steps, where one step corresponds 
to 0.0015 time units r. In the flow-aligned case, the sys- 
tem remained flow-aligned throughout the simulation. In 
the log-rolling case, two things happened during the ini- 
tial "equilibration" time: The director of the nematic 
phase turned very slowly into the flow direction, and the 
density of the nematic phase decreased slightly (implying 
that the thickness of the paranematic slab shrank) . After 
the initial "equilibration" time, the system had reached 
a state where half of the nematic slab was oriented at 45 
degrees to the flow, and the director of the other half was 
still perpendicular to the flow (see Fig. ^ . This state did 
not change any more over the next 2 million MD steps, in 
which the data were collected. Thus we could analyze the 
properties of this long-lived, yet presumably mctastable 
state and compare them with the properties of the stable 
flow-aligned interface. 

Since we use periodic boundary conditions, the inter- 
faces can move freely through the system. Moreover, 
they exhibit capillary wave fluctuations,i2iSliaS& which 
broaden the apparent profiles and makes their width de- 
pendent on the lateral system size. In order to eliminate 
these effects and determine local interfacial profiles, we 
adopt a block analysis technique introduced in our earlier 
studies of equilibrium interfaces^: We split the system 
of size L X Ly X L into columns of block size B X Ly X B. 
In the y direction, the columns are further divided into 
100 bins. Then we compute the local order tensor— in 
each bin. 
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FIG. 2: Director orientation profiles in the final configura- 
tion of the interface set up in tfie (a) fiow-afigned and (b) log- 
rolling geometry, with A*" = 115200 particles. In the isotropic 
slab, which is situated between y = — SOao and y = 0, the di- 
rection of the director is not well-defined. The symbols show 
the squared component of the director in different direc- 
tions: the shear gradient direction (plus, a = y), the director 
direction in the initial setup (closed circles), and the remain- 
ing direction (open diamonds). After the initial "equilibra- 
tion" time, these profiles stay roughly constant throughout 
the simulation. 



where I is the 3x3 unity matrix. The nematic order 
parameter S{y) in the bin centered at y is the highest 
eigenvalue of Q in that bin. From the profile S{y), we 
determine the local positions y^i and yjN of the two 
interfaces in the column under consideration, following 
the procedure described in Reference Then we calcu- 
late profiles for all quantities of interest and shift them 
by the amount y^j or yiN, respectively. This allows to 
perform averages over local profiles. Here, profiles from 
Nl-intcrfaccs are mirrored at y^i such that the nematic 
slab is always on the right side. The two interfaces from 
the log-rolling setup are discussed separately, since they 
were different (cf. Fig. The left interface, where the 
director on the nematic side is perpendicular to the flow, 
is labelled "log rolling" , and the right interface, where 
the director is oriented at 45 degrees with respect to the 
flow, is labelled "45 degrees". In the flow-aligned case. 
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we average over both interfaces. The interface positions 
y/jv and umi themselves are used to analyze the capillary 
wave spectrum of the interfaces. 

We should note that the two interfaces in the system 
are slightly correlated, due to the finite width of the 
isotropic and the nematic slab separating them. In par- 
ticular, the clastic interactions mediated through the ne- 
matic slab, which already affect the capillary wave spec- 
trum of a single interface will also cause long-range 
interactions between the interfaces. A detailed finite-size 
analysis, and data with very low statistical error, would 
be necessary to investigate these effects in detail. Unfor- 
tunately, this would require such an immense computa- 
tional effort, that it is not feasible for us at the moment. 

III. RESULTS 

We first present our results for the interface at the 
strain rate 7 = 0.001/t. Some color snapshots of the 
system can be found in a recent preliminary report for the 
John-von-Neumann institute in Jiilich, which is available 
online f§i 

The stress tensor, equal to the negative pressure tensor 
P, can be calculated via 

I I ^ iv \ 

cr = -P = -— ^ niiVi -f ^ ® % . (5) 

The shear stress is given by {(Txy + cTyx)/'2.. For the flow- 
aligned setup, we obtained 0.7 ± 0.1 • 10~"'eo/o'Q, and for 
the log-rolling setup, 1.3 ±0.1 • IQ-^eo/cTo. The fact that 
the shear stress is higher in the log-rolling setup {i.e., the 
average shear viscosity is higher) is another indication 
that the log-rolling state is metastable. The antisym- 
metric stress Gxy — (Jyx, as well as all other nondiagonal 
elements of cr, was zero within the error. For comparison 
with the equilibrium interface, we have also calculated 
the quantity 

l=^{Pyv-\iP^^+Pzz))- (6) 

In an equilibrium configuration with two interfaces in 
the {xz) plane, 7 is the interfacial tension. The equilib- 
rium value in our system was found to be^ 7 = 0.009 ± 
O.OOSeotTQ. Under shear, we obtain 7 = 0.003 ± 0.001 eo^o 
in the flow-aligned setup, and 7 = 0.015 ± 0.001 eoUo in 
the log-rolling setup. The difference from the equilibrium 
value can be attributed to two effects: First, the inter- 
facial properties change, and second, the diagonal com- 
ponents of the pressure in the bulk are no longer equal 
under shear. On principle, these two contributions can 
be separated, either by systematically varying the system 
size in the directions parallel and perpendicular to the in- 
terface, or via a detailed analysis of the profiles of 74i 
Unfortunately, the flrst method is computationally too 
expensive, and the statistics of the pressure profllcs was 
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FIG. 3: Local order parameter profiles at average shear rate 
7 = 0.001/r in a system with A'^ = 115200 particles for the 
flow aligned setup (solid), and the two interfaces from the 
log rolling setup (dashed and dot-dashed), compared to the 
corresponding equilibrium profile (dotted). Inset shows the 
difference AS between the profiles of the sheared and the 
equilibrium interface. 



not sufficient to obtain results from the second method. 
Therefore, we have not been able to separate the different 
contributions to 7, and the interpretation of our result is 
difficult. The difference between the pressure tensor com- 
ponents Pxx — Pzz is — 3± 1 • 10~^eoO'o in the flow-aligned 
setup, and zero within the error in the log-rolling setup. 

The local interfacial profiles have been analyzed fol- 
lowing the procedure described in the previous section, 
with the lateral block size B — 37.5cro (i. e., , the sys- 
tem was divided into 4 x 4 = 16 columns). Figs. |31 and 
0] show profiles of the order parameter and the density 
at the interface for the flow-aligned and the log-rolling 
setup, and compare them with the corresponding equilib- 
rium profiles. The order parameter profiles demonstrate 
that the interface broadens slightly under shear. In the 
log-rolling case, the effect is more pronounced than in the 
flow-aligned case. Moreover, the density of the coexisting 
nematic phase in the log-rolling case is slightly reduced. 

Next we examine the streaming velocity profiles at the 
interface, shown in Figs.ElandEl The streaming velocity 
profile in the flow direction clearly exhibits two different 
slopes: The local shear rate is smaller in the paranematic 
phase than in the nematic phase. Thus we observe shear 
banding. 

Even more remarkably, the shear also induces flow in 
the vorticity direction Vz (Fig. |SJ). The direction of that 
flow in the nematic and the paranematic phase is oppo- 
site. To our best knowledge, such a behavior has not 
been reported for homogeneous systems. (Even if it had 
been observed, it could be removed by a simple Galilei 
transformation.) Hence our effect is clearly induced by 
the interface. More precisely, the interface seems to in- 
duce a flow gradient in the vorticity direction. In the 
flow-aligned setup, where both interfaces are symmctri- 
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FIG. 4: Local density profiles at 7 = 0.001/r (iV = 115200 
particles). 



FIG. 5: Streaming velocity profiles in the flow direction, 
shifted such that = at y ^ 0, at ■y ^ 0.001/r (iV = 115200 
particles). The interface is located at y = 0. Profiles from NI- 
interfaces are mirrored at y = and Vx = 0. 



cal, the velocity gradient dies off far from the interfaces 
and the nematic and the isotropic slab move as blocks in 
opposite directions. In the log-rolhng setup, the effect is 
even stronger (by a factor of three), but the situation is 
complicated by the fact that the two interfaces are differ- 
ent, and the integrated effect of the flow gradient at the 
two interfaces is not compatible with the periodic bound- 
ary conditions. Therefore, we only show the data for the 
flow-aligned case here. 

In the direction of shear gradient (vy), no flow was 
observed. Likewise, the streaming angular velocity was 
indistinguishable from zero within the error (0.0002). In 
fact, Q cannot be strictly zero in the isotropic phase, 
since general considerations show that it should be one 
half the local shear ratei^ In our system, however, the 
local shear rate in the isotropic slab was 7 w 0.00025 
(cf. Fig. O, hence the resulting = 7/2 is within the 
statistical error. 

The fluctuations of the interface positions can be char- 
acterized by the capillary wave spectrum: We assume 
that we can parametrize the local position of an interface 
in the (x, z) plane by a single- valued function y = h{x, z). 
The inverse of the fluctuations of its two dimensional 
Fourier transform, (|/i(q)|^)~^, contains information on 
the local intcrfacial free energy. For example, an interface 
whose fluctuations are solely penalized by the interfacial 
tension 7 has the capillary wave spectrumS^ 

^bT ^ 2 

(IMq)P) " ^' 

The capillary waves of the equilibrium nematic/isotropic 
interface in our system have been investigated in detail in 
earlier workii^ The spectrum was found to be anisotropic, 
capillary waves being stronger in the direction perpendic- 
ular to the director than in that parallel to the director. 
This was recently explained theoretically by Elgcti and 
one of us within a Landau-dc Genncs treatment JS, Ac- 
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FIG. 6: Streaming velocity profile in the vorticity direction 
at 7 = 0.001/r for the flow-aligned setup. 



cording to the theory, the spectrum should be isotropic 
in the long-wavelength limit g — > 0. However, this limit 
is only reached on length scales of several thousand cor- 
relation lengths. On smaller length scales, the spectrum 
is dominated by strongly anisotropic cubic q'^-term and 
even higher order terms. 

The capillary wave spectrum is the result of a subtle 
interplay between the local director anchoring at the in- 
terface and the long-range elastic interactions. Since the 
latter are affected by shear, we expected that the spec- 
trum should change under shear. We have performed a 
capillary wave analysis as described in the previous sec- 
tion with a lateral column size B = IS.Sctq {i. e., the 
system was split into 8 x 8 = 64 columns). The re- 
sults are shown and compared with the spectrum of the 
equilibrium interface in Fig. Somewhat surprisingly. 
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FIG. 7: Capillary wave spectrum of the interface under 
shear with 7 = 0.001/r (iV = 115200 particles) in the flow- 
aligned case (a) and the log-rolling case (b), for wave vec- 
tors pointing in the vorticity direction and in the flow di- 
rection. Dashed lines and symbols show for comparison the 
equilibrium spectrum for wave vectors parallel to the direc- 
tor (squares, lower curve), and perpendicular to the director 
(circles, upper curve). 



shear does not have much noticeable effect on the capil- 
lary wave amplitudes. At low g-values, corresponding to 
large wave-lengths, the amplitudes (|/i(q)p) are higher 
than those of the equilibrium interface. This indicates 
that shear might reduce the effective interface tension, 
and is compatible with the observation that the inter- 
facial width broadens under shear. In the flow-aligned 
case, the capillary waves seem slightly enhanced in the 
vorticity direction. However, the effect is weak, and the 
deviations from the equilibrium values arc still within the 
error bars. 

Finally in this section, we discuss the phase coexis- 
tence and the stability of the interface at other strain 
rates. To study this, we have carried out a number of 
simulations of smaller, flow-aligned systems, where the 
length of the simulation box and the average density of 
particles was the same as before, but the lateral size was 
reduced by a factor of four (7200 particles in a box with 
aspect ratio 1:8:1). The interfaces remain stable up to 
an average strain rate of roughly 7 = 0.006/t. Fig. |S1 
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FIG. 8: Phase diagram of the flow-aligned system. The 
dashed "tie" lines connect coexisting states. I denotes the 
isotropic (paranematic) region, and N the nematic region. 
The data were collected with systems of size A'^ = 7200 parti- 
cles. 



shows the corresponding nonequilibrium phase diagram. 
Somewhat surprisingly, the densities of the two coexisting 
states are almost independent of the shear rate. This un- 
derlines our earlier statement that the phase separation 
in our system is almost exclusively driven by thermody- 
namic forces. The local shear rate is always higher in the 
nematic phase, i. e., the system is always shear thinning 
upon ordering. In contrast to predictions from theoret- 
ical model calculations the coexistence region does 
not close up. Rather, the interface disappears abruptly 
beyond the average shear rate 7 = 0.006/r. Coexistence 
regions similar to ours have been predicted for coexist- 
ing paranematic and log-rolling statesi^ Why we observe 
this structure here for paranematic states coexisting with 
flow-aligned nematic states is unclear, and will have to 
be the subject of future studies. 

Figs . 1^ and 1101 show the time evolution of interface de- 
struction at shear rates 7 = 0.1/t and 7 ~ 0.01/r. The 
initial configuration is an equilibrium interface. At the 
shear rate 7 = 0.1/t, the shear aligns the liquid crystal 
throughout the paranematic slab, thereby destroying the 
difference between the two phases. This is the primary 
process and happens on a time scale of roughly 30r. In a 
second step, the density slowly adjusts. The density pro- 
file therefore persists much longer than the order param- 
eter profile (Fig. |51b). At the shear rate 7 = 0.01/t, the 
destruction process is less homogeneous and nucleates at 
the interfaces: They slowly move inwards and gradually 
lose their structure (Fig. I10|l . This happens on a time 
scale of roughly 300r. The time scale of destruction at 
average shear rate 0.01/r is thus about ten times as large 
as at shear rate 0.1/r. We conclude that the time scale of 
destruction is roughly proportional to the time scale 1 /-j 
introduced by the shear. This underlines the necessity 
of performing very long simulation runs at lower shear 
rates. In particular, it is important to note that at the 
shear rate 7 = 0.001/r, where we have performed most 
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FIG. 9: Destruction of the interface in the flow ahgned system 
at 7 = O.l/r (A^ = 115200 particles), a) Order parameter 
profile b) Density profile. 



of our simulations, the interface showed no sign of desta- 
bilization after a run of length 6000t. Therefore we can 
deduce that the interface is indeed stable. 



IV. CONCLUSIONS AND OUTLOOK 

We have presented an extensive nonequiUbrium molec- 
ular dynamics simulation of liquid crystal interfaces in 
shear flow, and compared the properties of these inter- 
faces with the structure of the corresponding equilibrium 
interface. Two situations were considered: In the first, 
flow-aligned setup, the interface was set up such that the 
director in the nematic phase pointed in the direction 
of shear gradient; in the second, log-rolling setup, the 
director pointed in the direction of the vorticity. The 
log-rolling state appeared to be unstable with respect to 
the flow-aligned state, but extremely long-lived, so that 
it could be analyzed as well. In both the flow-aligned 
and the log-rolling system, the two coexisting states cx- 



FIG. 10: Destruction of the interface in the flow aligned 
system at 7 = 0.01/r (TV = 115200 particles). 



perience the same, common, stress. The third possible 
situation, phase coexistence at common strain (Fig. ^ 
right), was not yet studied. This will be done in future 
work. 

We found that the structural properties of the interface 
are not dramatically affected by the presence of shear. 
The interfacial width increases slightly, and the capil- 
lary wave amplitudes are slightly enhanced at large wave- 
lengths. No other noticeable effects were observed. For 
example, the biaxiality was still basically zero (data not 
shown) . 

In contrast, the streaming velocity profllcs revealed a 
number of interesting features: First, we observe shear 
banding, i. e., the shear profile is inhomogcncous. The 
local shear rate in the nematic phase is distinctly higher 
than in the paranematic phase. Second and surprisingly, 
the interface also induces a streaming velocity gradient 
in the vorticity direction. As a result, vorticity fiow in 
opposite directions is induced in the paranematic and 
the nematic phase. This flow is symmetry breaking, it 
destroys the mirror symmetry at the (2:y)-plane. The 
effect is small but significant. So far, we have no expla- 
nation. It has not been predicted by theory; however, the 
theoretical studies of shear banding cited in the introduc- 
tion were all set up such that the vorticity fiow was zero 
by construction. Wc arc not aware of any fundamental 
argument against nonzero vorticity flow. Hence our re- 
sult is perplexing, but not impossible. We hope that it 
will stimulate further theoretical work on shear banding, 
and that the mechanism and the conditions for symmetry 
breaking will be clarified. 

The effect is not only of academic interest, it could also 
be useful from a practical point of view. For example, if 
one component enriches in one phase, and the other in 
the other phase, shear could possibly be used to separate 
the two components. 

The present study has focused on interfaces which are 
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mainly stabilized by thermodynamic driving forces. As 
mentioned in the introduction, nonequilibrium interfaces 
can also be stabilized mechanically. In that case, the in- 
tcrfacial properties presumably differ much more strongly 
from those of equilibrium interfaces than in the present 
case. We hope that we will succeed in stabilizing and 
studying such an interface in the future. 



N = 7200 systems within an undergraduate physical 
chemistry project. The parallel MD program GBmcga 
used in this work was originally developed by the EP- 
SRC Complex Fluid Consortium. This work was funded 
by the German Science Foundation, and the simulations 
were carried out at the John-von-Neumann computing 
center in Jiilich. 
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